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just above the threshold. We also consider the bound-state pole condition in finite volume, based 
on Liischer's phase-shift formula and then find that the condition is fulfilled only in the infinite 
volume limit, but its modification by finite size corrections is exponentially suppressed by the 
spatial lattice size L. These theoretical considerations are also numerically checked through lattice 
simulations to calculate the positronium spectrum in compact scalar QED, where the short-range 
interaction between an electron and a positron is realized in the Higgs phase. 
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1. Introduction 

Signatures of bound-state formation in finite volume are of main interest in this paper. In the 
infinite volume, the bound state is well defined since there is no continuum state below threshold. 
However, in a finite box on the lattice, all states have discrete energies. Even worse, the lowest 
energy level of the elastic scattering state appears below threshold in the case if an interaction is 
attractive between two particles Therefore, there is an ambiguity to distinguish between the 
shallow (near-threshold) bound state and the lowest scattering state in finite volume in this sense. 

We may begin with a naive question: what is the legitimate definition of the shallow bound 
state in the quantum mechanics? In the scattering theory, poles of the ^-matrix or the scattering 
amplitude correspond to bound states iQ]. It is also known that the appearance of the 5- wave bound 
state is accompanied by an abrupt sign change of the S-wave scattering length It is interpreted 
that formation of one bound-state raises the phase shift at threshold by n. This particular feature 
is generalized as Levinson's theorem [^. Thus, it is interesting to consider how the formation 
condition of bound states is implemented in Liischer's finite size method, which is proposed as a 
general method for computing low-energy scattering phases of two particles in finite volume [jl|]. 

In this paper, we discuss bound-state formation on the basis of the Liischer's phase-shift for- 
mula and then present our proposal for identifying the shallow bound state in finite volume. To 
exhibit the validity and efficiency of our proposal, we perform numerical studies of the positron- 
ium spectroscopy in compact scalar QED model. In the Higgs phase ofU{\) gauge dynamics, the 
photon is massive. Then, massive photons give rise to the short-ranged interparticle force between 
an electron and a positron, which is exponentially damped. In this model, we can control positro- 
nium formation in variation with the strength of the interparticle force and then explore distinctive 
signatures of the bound-state formation in finite volume. The contents of this paper are based on 
our published work 

2. Bound-state formation in Luscher's formula 

In quantum scattering theory, the formation condition of bound states is implemented as a 
pole in the S-matrix or scattering amplitude. Here, an important question naturally arises as to how 
bound-state formation is studied through Luscher's phase-shift formula Intuitively, the pole 
condition of the ^-matrix: S = e^'*'^' = c°||||p|^) is expressed as 

cot5Q{p) = i, (2.1) 

which is satisfied at /j^ = — where positive real y represents the binding momentum. In fact, as 
we will discuss in the following, such a condition is fulfilled only in the infinite volume. However 
the finite-volume corrections on this pole condition are exponentially suppressed by the size of 
spatial extent L. 

It was shown by Liischer that the 5-wave phase shift 5q can be calculated by measuring the 
relative momentum of two particles p in a. finite box L} with a spatial size L through the relation 

71:3/2^^ 

tango(p) = ^ OS atq = Lp„/27i, (2.2) 
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where the generalized zeta function ^oo{s,q^) = 



An ■ 



-neZ 



an 



-q^) ^ is defined through analytic 



continuation in s from the region 5^>3/2tos = l [|T|. For negative q^, an exponentially convergent 
expression of the zeta function ^^^{s^q^) has been derived in Ref. iQ]. For s = 1, it is given by 



iroo(l,^^) = -7rV2yZ^+£'^, 



(2.3) 



where Ln^z^ ni^ans the summation without n = (0,0,0). We now insert Eq. {13) into Eq. {12) and 
then obtain the following formula, which is mathematically equivalent to Eq. (2.2) for negative q^: 



neZ' 



(2.4) 



The second term in the r.h.s. of Eq. (2.4) vanishes in the limit of q^ — oo. It clearly indicates 
that negative infinite q^ is responsible for the bound-state formation. Therefore, in this limit, the 
relative momentum squared approaches —7^, which must be non-zero. Meanwhile, the negative 
infinite q^ turns out to be the infinite volume limit. 

According to the original paper [|l|], for negative q^, we introduce the phase Oq{k), which is 
defined by an analytic continuation of 5q into the complex p plane through the relation tanao(K") = 
— /tan5o(/j), where k = —ip. As a result, the bound-state pole condition in the infinite volume 
reads cot 09(7) = — 1 for the binding momentum 7 ||l|]. Then, Eq. ( ^i^l ) can be rewritten in terms of 
the phase ao as 



lim cotaofjc) 



,-s/vLy . 



v=i VvLy 



-1 + 



Lj 



e-^^+. 



(2.5) 



where the factor A^v is the number of integer vectors n S with v = n^. Therefore, it is found that 
although a bound-state pole condition is fulfilled only in the infinite volume limit, its modification 
by finite size corrections is exponentially suppressed by the spatial extent L in a finite box \ We 



can learn from Eq. (2.5) that "shallow bound states" are supposed to receive larger finite volume 
corrections than those of "tightly bound states" since the expansion parameter is scaled by the 
binding momentum 7. 



3. Novel view from Levinson's theorem 



If the 5-wave scattering length ao, which is defined through ao = liiUp^o tan 5o (/?)//?, is suffi- 
ciently smaller than the spatial size L, one can make a Taylor expansion of the phase-shift formula 



(2.2) around q^ = 0, and then obtain the asymptotic solution of Eq. (12). Under the condition 
<C where pL represents the reduced mass of two particles, the solution is given by 



27rao 



do 

l+Cl-^+C2 



+ 



(3.1) 



which corresponds to the energy shift of the lowest scattering state from the threshold energy. The 
coefficients are c\ = —2.837297 and C2 = 6.375183 [II]]. An important message is received from 



^Although it was pointed out how the bound-state pole condition could be implemented in his phase-shift formula 
in the original paper [j^], this important fact has been firstly reported in Ref. |^. 
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Eq. (3.1). The lowest energy level of the elastic scattering state appears below threshold {AE < 0) 
on the lattice if an interaction is weakly attractive (ao > 0) between two particles. This point makes 
it difficult to distinguish between near-threshold bound states and scattering states on the lattice. 

Here, a crucial question arises: once the S-wave bound states are formed, what is the fate of 
the lowest S-wave scattering state? The answer to this question might provide a hint to resolve 
our main issue of how to distinguish between "shallow bound states" and scattering states. A 
naive expectation from Levinson's theorem in quantum mechanics is that the energy shift relative 
to a threshold turns out to be opposite in comparison to the case where there is no bound state. 
Levinson's theorem relates the elastic scattering phase shift d[ for the Z-th partial wave at zero 
relative momentum to the total number of bound states (Ni) in a beautiful relation 5/(0) = NiH \ 
Therefore, if an 5-wave bound state is formed in a given channel, the S'-wave scattering phase 
shift should always be positive at low energies. This positiveness of the scattering phase shift 
is consistent with a consequence of the attractive interaction. Conversely, the 5-wave scattering 
length may become negative (ao < 0), especially for the shallow bound-states [Q]. Consequently, 
according to Eq. (p]), possible negativeness of the scattering length gives rise to a positive energy- 
shift of the lowest scattering state relative to the threshold energy. In other words, the lowest 
scattering state is pulled up into the region above threshold. Therefore, the spectra of the scattering 
states quite resembles the one in the case of the repulsive interaction. If it were true, we can observe 
a significant difference in spectra above the threshold between the two systems: one has at least 
one bound state (bound system) and the other has no bound state (unbound system). 



4. Numerical results 



To explore signatures of bound-state formation on the lattice, we consider a bound state 
(positronium) between an electron and a positron in the compact QED with scalar matter : 

5sQED[t/,0,»F] =/3 £ [1 -9?{[/,,^v}] -hY,^{^lU,^^^,+^] + Y,^MU\,^y^y, (4.1) 

plaq. link sites 

where jS = 1/e^ and the constraint \(^x\ = 1 is imposed. This action is described by the compact 
U{\) gauge theory coupled to both scalar matter (Higgs) fields and fermion (electron) fields *P. 
In this study, we treat the fermion fields in the quenched approximation. 

Our purpose is to study the S'-wave bound state and scattering states through Liischer's finite 
size method, which is only applied to the short-ranged interaction case. Thus, we fix j8 = 2.0 and 
/j = 0.6 for the compact ?7(1)-Higgs action to simulate the Higgs phase ofU{\) gauge dynamics, 
where massive photons give rise to the short-ranged interparticle force between an electron and 
a positron. We generate U{\) gauge configurations with a parameter set, (j8,/z) = (2.0,0.6), on 
X 32 lattices with several spatial sizes, L = 12, 16,20,24,28 and 32. Details of our simulations 
are found in Ref. [|]. 

Once the parameters of the compact ?7(1)-Higgs part, (j8,/z), are fixed, the strength of an 
interparticle force between electrons should be frozen on given gauge configurations. However, 
if we consider the fictitious Q-charged electron, the interparticle force can be controlled by this 

^ Strictly speaking, this form is only valid unless zero-energy resonances exist. 



4 



Identification of shallow two-body bound states infinite volume 



Shoichi Sasaki 



> ground 

■ 1 St excited 

► 2nd excited 



Q=3 L=28 



> ground 

1 1 St excited 

► 2nd excited 



I 



Q=4 L=28 



Figure 1: The effective mass plots for each eigenvalue of the transfer matrix in the '5o channel on the lattice 
with L — 28. Full circles, squares and diamonds represent the ground state, the first excited state and the 
second excited state. The left (right) panel is ior Q = 3 (Q — 4). 



charge Q since the interparticle force is proportional to (charge Q)^. Within the quenched ap- 
proximation, this trick of the 2-charged electron is easily implemented by replacing U{1) link 
fields as Ux^^ — > U^^ = Hf^^Ux^^ into the Wilson-Dirac matrix According to our previous 
pilot study 10], numerical simulations are performed with two parameter sets for fermion (electron) 
fields, {Q,k)={3, 0.1639) and (4, 0.2222). As we will see later, the former case (Q = 3) corre- 
sponds to the unbound system, while the latter case {Q = 4) corresponds to the bound system where 
the positronium state can be formed. Here K, which is the hopping parameter of the Wilson-Dirac 
matrix, is adjusted to yield almost the same electron masses Me ~ 0.5 for both charges. 

We are especially interested in the '^o and ^Si states of the system, where the electron- 
positron bound state (positronium) could be formed even in the Higgs phase. ^Sq and ^Si positro- 
nium are described by the bilinear pseudo-scalar operator ^xYs^x and vector operator ^xY^i^x 
respectively. Therefore, we may construct the four-point functions of electron-positron states 
based on the above operators. We are interested in not only the lowest level of two-particle 
spectra, but also the 2nd and 3rd lowest levels. In order to extract a few low-lying energy lev- 
els of two-particle system, we utilize the diagonalization method [^. We consider three types 
of operators for this purpose: Q.p{t) = L-^l^^'¥{x,t)r^{x,t), Q.wit) = L-^l^^yW{y,t)r^'{x,t) 
and auit) = L-^^^^yW{y,t)r¥{x,t)e'P'<''-^'^ where pi = ^(1,0,0) and T = 75 (7^) for the '^o 
(^^i) e^e^ state. We construct the 3x3 matrix correlator from above three operators Gij{t) = 
(0|n,(f)nj(0)|0) and then employ a diagonalization of a transfer matrix. As shown in Fig. [T], the 
diagonalization method with our chosen three operators successfully separates the first excited state 
and the second excited state from the ground state [Q]. 

4.1 Sign of energy shift 

In Figs. |2[ we show energies of the ground state and also excited states in the system as a 
function of spatial lattice size L. The dashed lines and curves represent the threshold energies 2Mg 
and 2Ee{p\), which are evaluated by measured energies of the single electron with zero momentum 
Po = ^(0,0,0) and nonzero lowest momenta pi = ^(1,0,0) respectively. Two left panels are for 
the ^5o channels, while the right panel is for the channel. 
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Figure 2: Energies of the ground state and excited states in the '5o (the left and middle panels) and ^Si (the 
right panel) channels of the e^e+ system as functions of spatial lattice size. The left figure is for Q = 3, 
while the middle and right panels are for Q = 4. 



Let us focus on results in the channel. The energy level of the ground state for 2 = 3 in 
the left panel appears close to the threshold. An upward tendency of the L-dependence toward the 
threshold energy is observed as spatial size L increases. This is consistent with a behavior of the 
lowest scattering state predicted by Eq. (3J_) for the weakly attractive interaction without bound 
states. On the other hand, we clearly see the presence of a bound state for Q = 4, which certainly 
remains finite energy gap from the threshold even in the infinite-volume limit. The most striking 
feature is our observed L-dependence of the energy level of the second lowest state for {Q = 4). 
Clearly, this energy level approach the threshold energy from above. The energy shift vanishes as 
the spatial size L increases. Therefore, the second lowest energy state must be the lowest scattering 
state with the repulsivelike scattering length {qq < 0). In addition, the level of the second lowest 
scattering state are located near and below (above) the threshold energy 2Eg{pi) for Q = 3 (Q = 4). 

In the channel for 2 = 4, the binding energy B is rather large as B Mg/2. The observed 
bound state should be a "tightly bound state" rather than a "shallow bound state". On the other 
hand, we observe that the bound state in the channel (the right panel) is much near the threshold 
energy. Although the ground state lies too close to the threshold energy to be assured of bound- 
state formation, the distinctive signature of bound state is given by an information of the excited 
state spectra. The second lowest state appears just above the first threshold 2Mg, but far from the 
second threshold 2£'e(pi). Therefore, we can conclude: the ^Si ground state should be the shallow 
bound state, of which formation clearly induces the sign of the scattering length to change [^]. 



4.2 Bound-state pole condition 

A rigorous way to test for bound-state formation would be to use an asymptotic formula for 
finite volume correction to the pole condition as Eq. (2.5). In Figs.^, we plot the value of cotao 
versus the spatial lattice extent L for either (left) and (right) channels. Full circles are 
measured value at five different lattice volumes. At first glance, we observe that the phase cot Gq 
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gradually approaches —1 as spatial lattice extent L increases for either channels. 



We next examine the L-dependence of cot Gq by reference to Eq. (2.5), where the finite volume 
corrections on the bound-state pole condition aie theoretically predicted. The solid and dashed 
curves represent fit results with a single leading exponential term and three (six) exponential terms 
in the ^Sq CSi) channel. All five data points are used for those fits in the ^Sq channel, while the 
four data points in the region 20 < L < 32 are used in the channel. The fitting with the three 
(six) exponential terms yields a convergent result of /in the ^Sq CS\) channel. Either fit curves in 
Figs. H reproduce all data points except for data at the smallest L in the ^Si channel. Therefore, we 
confirm that the ground state in the channel at least for L > 20 can be identified as a shallow 
bound state without ambiguity. 

5. Summary and conclusion 

In this paper, we have discussed formation of an S-wave bound-state in finite volume on the 
basis of Liischer's phase-shift formula. We have first showed that although a bound-state pole 
condition is fulfilled only in the infinite volume limit, its modification by the finite size corrections 
is exponentially suppressed by the spatial extent L in a finite box L^. We have also confirmed that 
the appearance of the S-wave bound state is accompanied by an abrupt sign change of the S-wave 
scattering length even in finite volume through numerical simulations. This distinctive behavior 
may help us to discriminate the shallow bound state from the lowest energy level of the scattering 
state in finite volume simulations. 
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